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ABSTRACT 

Although the magnetorotational instability (MRI) has been widely accepted as a pow- 
erful accretion mechanism in magnetized accretion disks, it has not been realized in the 
laboratory. The possibility of studying MRI in a rotating liquid-metal annulus (Couette 
flow) is explored by local and global stability analysis. Stability diagrams are drawn 
in dimensionless parameters, and also in terms of the angular velocities at the inner 
and outer cylinders. It is shown that MRI can be triggered in a moderately rapidly 
rotating table-top apparatus, using easy-to-handle metals such as gallium. Practical 
issues of this proposed experiment are discussed. 
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1 INTRODUCTION 

Astrophysical magnetic fields have long been recognized to 
be important but difficult to un derstand. A prominent ex- 
ample concerns accretion di sks ( Bhakura & Sunyaev 197J ; 



Lynden-Bell & Pringle 1974), where orbiting plasmas grad- 



instability of a magnetized but Rayleigh-stable fiuid, the 
mag netorotational i nstabilitv (MRI), was discovered decades 
ago (Velikhov 195!:; Chandrasckhar 196C) but did not come 
to the attention of the astrophyisical commu nity until re- 
cently rediscovered (Balbus & Hawley 1991a), despite gen- 



ually accrete onto a central mass. Three types of central 
object occur: protostars (stars in formation), collapsed stars 
in binary systems (white dwarfs, neutron stars, and black 
holes), and supermassive black holes in active galactic nuclei 
(quasars et al.). In addition to the accretion process, jets and 



eral reco gnition that magnetic effec ts might somehow be im- 
portant (Shakura & Sunyaev 197J). Since then, many ana- 
lytic and numerical studies of the MRI have been performed 
under increasingly complex and realistic assumptions, in- 
cluding such effects as finite resistivity, global boundary con- 



othHi spe cLacLilai plienuincna may be magne lically diivenby 
the 
the 



ditions, and nonlinearity in two and three dimensions (Bal 



diak (e.g. Meier, Koidc, fc Uchida 2001). Undcratanding 



dyiiamicb and dissipaLiun mcchaiiibiub in accicLiuii disks 



bus fc Hawley 1991b| ; |Curry, Pudritz, fc Sutherland 1994 
Blaes fc Balbus 1994|; Brandenburg et al. 199E|; |Matsumotq 



holdj au impuitauL kL_y Lu uiidLuUiidiug mmy iicLiuc iuLiu" 



^^imal995; Hawley, Gammie, fc Balbus 1996 ; Stone et 



nor hical avatcma 



yholc. 



The accretion rates cannot be due to ordinary molecu- 
lar or plasma viscosity because of the extraordinarily high 
Reynolds numbers involved. Theorists hav e often appealed 
to hydrodynamic turbulence (Pringle 1981), but recent nu- 



al. 1996; Gammie 1996; lin 199£; Sano fc Miyama 1999 



Fleming, Stone, fc Hawley 2000^ [Hawley 2000^ 

Despite its popularity and importance, however, the 
MRI has never been realized in the laboratory or demon- 
strated observationally. Laboratory plasma experiments are 
primarily magnetically driven, and the observed fiows, often 
induced as secondary effects of other instabilities, are small 
compared to the Alfven speed. On the other hand, the exist- 
ing body of experim ental work on magnetized Couette flow 
using liquid metals (Donnelly fc Ozima 196C ; Donnelly fc 



merical simulations indicate that nonmagnetic disks are sta- 



bilized by their positive ang ular momentu m gradient ( Bal 
bus, Hawley, fc Stone 1996| |Cabot 199^ ; [Hawley, Balbus 



^ Winters 199i: ); in effect, Rayleigh's stability criterion ap- 
pea fa to aupprcaa loo al nonaxiaymmc trio aa well aa axiaym 



mei ric disturbances (Rayleigh 1916) . Linear axisymmetric 



Ozima 1962 



Donnelly fc Caldwell 1964| ; [Brahmc 1970D has 
focused on magnetic stabilization of the Ray leigh instabil- 
ity, as first analyzed by Chandrasckhar (1961). In this Let- 
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ter, we explore the feasibility of a Couette-flow experiment 
dedicated to MRI. 
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2 LOCAL STABILITY ANALYSIS 



tions lead to the following dispersion relation: 



Couette flow inv olves a liquid confined between rotating 
coaxial cylinders (Couette 189C). Let their radii be r\ < r2, 
and their angular velocities 572- In steady state, the radial 
angular momentum flux, h-2iTr- pv -r^ {—dO,/ dr), is constant 
with radius, where h is the depth of the liquid, p is its den- 
sity, and V is its kinematic viscosity. If dh/dr — 0, then the 
angular velocity of the liquid satisfies r^dQ/dr =const., so 
that 



n(r) 



(1) 



where a = (ri2r-i — Qirl)/{r2 — rf) and b = rir2{^li 



^2)/{r2 — rl). The Rayleigh stability criterion is aQ > 0. 

The dynamics of liquid metals is well described by 
the incompressible and dissipative magnetohydrodynamic 
(MHD) equations. 
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where V is velocity, B is magnetic field, is magnetic dif- 
fusivity, and p is a scalar potential incorporating both pres- 
sure and gravity. In cylindrical coordinates, the equilibrium 
quantities are Bq — (0,0, B) and Vq — {0,rQ,,0), and the 
balance of forces is dpo/dz = and dpo/dr — prO? . 

WKBJ methods describe the stability of this system 
very well even on the largest scales. Using cylindrical co- 
ordinates, the perturbations are B\ = {Br, Bg, B^) and 
Vi — {Vr, Vg, Vz), all proportional to exp(7t — ik^z — ikrv), 
so that 7 is the growth rate and the perturbations are ax- 
isymmetric. The minimum fcz and kr are assumed to be 7r//i 
and 7r/(r2 — ri), respectively, so that the total wavenum- 
ber k = y/k'i + kr — fczVl + e^, where e = h/{r2 — ri) is 
the elongation of a toroidal cross-section of the liquid metal 
annulus. The linearized equations of motion are 
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where the epicyclic frequency is defined by k? = 
{l/r^)d{r*0'^)/dr = 4^^ + dO^ /d\nr and pi is the per- 
turbed pressure. The vertical induction equation is not 
needed since Bz can be deduced from V-B = 0. These equa- 
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The Alfven speed is Va = B/^/pop. This dispersion rela- 
tion is identical to the one derived for accretion disks in the 



incompressible limit (Sano & Miyama 1999). 

Introducing a dimensionless vorticity parameter, C, = 
{l/rQ)d(r^n)/dr ^ 2 + d\nQ,/d\nr, we have = so 
that the Rayleigh stability criterion becomes C ^ 0- Simi- 
larly, dO'^/dlnr = 2n^(C - 2). There are three other rele- 
vant frequencies: resistive, cj^ = "qk^; viscous, ui^, = vk"^; and 
Alfvenic uja = \kzVA\- Because liquid metals are far more 
resistive than viscous, ujr) serves as a base frequency in the 
following three dimensionless parameters: magnetic Prandtl 
number, Pm = w^/tJ,,; Lundquist number, S = LOA/i^r]\ and 
magnetic Reynolds number, Ria = 0/uir,- The astrophysical 
literature gives several inequivalent definitions of "magnetic 
Reynolds number," some corresponding to our S. Some in- 
volve the sound speed, Cs, for although MRI is essentially 
noncompressive, vertical force balance in an astrophysical 
disk relates Cs to the half-thickness: Cs ~ hO. The free en- 
ergy for MRI derives from differential rotation, represented 
in our dimensionless system by a combination of ( and Rm', 
but magnetic field, represented by S, is required to liberate 
this energy from hydrodynamic constraints. 

Using the normalized growth rate, 7/0;,, —> 7, the dis- 
persion relation can be rewritten as 

[(7 + ^'m)(7 + l) + S^]'(l + f') 

+2(:rI{j + if -2(2- ORlS^ = 0. 

It can be shown that there are no purely imaginary roots for 
7 as follows. Suppose that 7 = ioj for real and nonzero (jj. 
From the imaginary part of the equation above, one finds 



S'+P„, + (l + Pn.)Ca 



Pn 



1 + Pn, 

Substituting this into the real part yields 



1 + 



This is a quadratic equation in ( which must be real. But 
the discriminant 



Pm 



is negative unless S = or Pm = 0. Hence there are no 
purely oscillatory modes. On the other hand, all roots of the 
dispersion relation for 7 have negative real parts as Pm 
(nonrotating fiow). Hence the transition to instability occurs 
through 7 = 0, and the necessary and sufficient condition for 
stability is that the value of the dispersion relation at 7 = 
remain positive as Pm increases: 



(Pm + SY(1 + e') -I- 2CPm - 2(2 - ORlS^ > 0, 



(2) 



which can be taken into various limits. 

Nonmagnetic limit. As r; — + 00, the three terms 5*, 
RmS^, and R^S^ approach zero faster than the others, lead- 
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Figure 1. Stability of a rotating liquid metal annulus in dimen- 
sionless parameter space of (5, Q) at a few values of i?m for the 
case of e = 1. Areas to right of the curves indicate stability. 



ing to the stability condition P^{1 + e^) + 2CRL > 0. Sta- 
bility occurs when ^ > 0, and also when < if the Taylo r 
number-like expression -2(n'^/u^k'^ < 1 + ( [Taylor 1923| ). 

Ideal MHD limit. As 0, the other two terms dom- 
inate, with stability for ^^(l + e^) > 2(2 - C)-Rm- Instability 



occurs at sufficie ntly weak fields (small S) unless C > 2, i.e. 
dlnD,/dlnr > ( [Balbus fc Hawley IQQl^ )- 

Small Pm limit. In liquid metals, usually viscosity is 



much smaller than resistivity, 
reduces to 

A.* 



10-^ AsPn, 



2S^ 



S2 + 1 2Rl{S^ + 1) 



0, eq. (H) 



(3) 
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Figure 2. Stability of a rotating liquid metal annulus in i?ni 
and space. Here the stability can be divided into 3 regions: 
region (I) is hydrodynamically unstable but can be stabilized by 
a finite magnetic field, as exemplified by point B. Region (II) is 
hydrodynamically stable but can be destabilized by presence of a 
magnetic field (MRI), as exemplified by points A and C. Region 
(III) is always stable. Results from global eigenmode analysis are 
also shown: dotted lines for conducting boundary conditions and 
dashed lines for insulating boundary conditions. 



Table 1. Parameters for a gallium annulus with r\ 
r2 = 0.15m, and h = 0.1m. 



0.05m, 



point 
A 
B 
C 



Qi (rpm) 
3600.00 
3600.00 
5089.77 



f22(rpm) 
435.00 
390.00 
620.70 



0.3319 
0.3143 
0.4715 



C 

0.06293 
-0.01899 
0.06984 



3 STABILITY DIAGRAMS AND GROWTH 
RATES 

The stability condition (^ defines a two-dimensional surface 
in the parameter space (5, Rm) at fixed e. To illustrate the 
dependence on these parameters, we vary only two of them 
at a time. Stability boundaries in the {S, Q plane are shown 
in Fig. 1 for the case of e = 1. When C < 0, the annulus is 
unstable hydrodynam ically to the Rayleigh mode at S = 
but can be stabilized ( Shandrasekhar 1961) by a large mag- 
netic field (large 5"). When > 0, the annulus is stable at 
zero field but unstable at some 5* > if Rm is large enough. 
Stability returns at even larger S. The unstable region ex- 
tends to larger 5* and at larger R n^. Stability at 5* = and 
as S ^ ex are hallmark s of MRI (Balbus & Hawley 1991a; 



Balbus & Hawley 1998). (In ideal MHD, instability extends 



formally to S = 0^.) It can be seen that there is a maximum 
^ above which MRI is absent for a given Rm- From eq.(ph. 



at 



the S value given by = y/l + 4Rl/{l + e'^) 



(4) 



- 1. 



Figure 2 shows stability in the (i?m, C) plane. Region (I) 
is hydrodynamically unstable but can be stabilized by a finite 
magnetic field. This region has be en extensively studied both 
theoretically and experimentally ( Chandrasekhar 1961), and 
is exemplified by point B. Region (II) is hydrodynamically 
stable but destabilized by a magnetic field. This is the MRI 
region, and has never been studied experimentally. Growth 
rates at points A and C are given below. Region (III) is 
always stable. The boundary between regions (II) and (III) 
is given by eq. (^. 

It is useful to project the stability diagram onto ex- 
perimentally controllable parameters. To apply the local 
dispersion relation, we take ^ = 2a/Q [from eq.(|l|)] and 
Rui = il/ujri, with fl = \^niQ2. Figure 3(a) shows stability in 
the plane (57i,572) for an annulus of dimensions ri=0.05m, 
7-2=0. 15m, and h—Q.lm (hence e = 1) filled with gallium 
(p ~ 6 X 10^kg/m^ 7? ~ 0.2mVs, u ~ 3x 10~^mVs)- Ta- 
ble 1 lists the physical parameters at points A, B, and C. 
The corresponding growth rates are shown as functions of 
magnetic field in Fig. 3(b). 

The applicability of WKBJ is subject to doubt, since the 
most unstable wavelengths are larger than the gap width and 
cylinder height. In fact, global analysis shows that the eigen- 
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Figure 3. Stability diagram of a rotating gallium annulus in 
Oi and 02 space with dimensions ri=0.05m, r2=0.15m, and 
h =0.1m. The growth rates of points A, B, and C, corresponding 
to those in Fig. 2, are also shown as functions of magnetic field in 
(b). Results from global eigenmode analysis are also shown: dot- 
ted lines for conducting boundary conditions and dashed lines for 
insulating boundary conditions. 



functions are nonsinusoidal and sensitive to the boundary 
conditions. Yet the growth rates are remarkably robust. A 
linearized, finite-difference, initial-value code was written to 
to detect the fastest growing mode. Periodic boundary con- 
ditions were used in z. Radial boundaries were impenetrable 
and no-slip {5Vi = 0), and electrically either perfectly in- 
sulating or perfectly conducting. Results are compared with 
the WKBJ analysis in Figs. 2 and 3. Figure 4 shows eigen- 
modes for the parameters of point C with 5=0.3 Tesla. Dif- 
ferences between the conducting and insulating cases can be 




Figure 4. Eigenmodes for conditions given by point C in Fig. 3 
at _B = 0.3 Tesla with conducting (a) and insulating (b) radial 
boundaries. Here, solid (dotted) lines represent positive (nega- 
tive) values; x ^i^d ij) are poloidal flux and stream functions, re- 
spectively. 



seen near the inner boundaries. A Hartmann layer (Hart 



within a radial thickness of ~ -^/vrj/VA < 1mm (not visible 
in Fig. 4), forms at the inner conducting boundary as the 
Lorentz force, 5jr x , balances with the viscous force. Nev- 
ertheless, the growth rates are remarkably similar to those 
of the local analysis, which therefore should suffice for pre- 
liminary experimental desig n. [Details of the glo bal analysis 
will be reported elsewhere (Goodman & Ji 2001).] 

A fully nonlinear incompressible MHD code has been 
developed to study the problem in three dimensions 



mam 1937), consisting of large toroidal and axial velocities (Kageyama, Ji, & Goodman 2001). Initial results from lin 
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ear and two-dimensional runs of this code with conducting, 
fre ely-slipping boundary conditions agree with the local and 



The experiment will be far more r esistive than most ac 



cretion disks, thou gh perhaps not all (Gammie 1996 



global analyses. For example, for the conditions given by the mie fc Menou 1998). Simulations indicate that when the field 



Gam- 



point C in Figs. 2 and 3 at _B = O.STesla, the growth rate is 
21.67 s~^ from simulations, 21.90 s~^ from global analysis, 
and 19.10 s^^ from local analysis. 



nT.SCTTSSTONS AND CONCTJISTONS 



is generated by the disk itself (magnetic dynamo), then the 
large-scale field is nearly toroidal, the important instabilities 
are nonaxisymmetric, and the turbulence sustains itself only 
if S and iim are much larger than the experiment proposed 
here will achieve { ^i^aus^^i^^^w\c^^^^\^QVo\ Brandcnburg_eq 
al. 1995 ; Hawley , Gammic, fc Balbus ano, Inutsuku 



SeVeral Yuau^ti muHt be explored before commiLLliig to aii ex- 
perimental design. The first is geometric optimization with 
regard to aspect ratio [A = (r2 + ri)/(r2 — ri)] and elon- 
gation [e = hjiji — ri)]. Obviously, e ~ 1 is desirable to 
minimize volume, and therefore expense, at a given growth 
rate. Less obviously, aspect ratios close to unity cause the 
eigenmodes and growth rates to be dominated by the in- 
ner cylinder, which is undesirable if the aim is to imitate 
relatively uniform local conditions within an astrophysical 
disk. Therefore, moderate values of 2) and e(~ 1) are 
preferred. 

The periodic vertical boundary conditions used in all of 
our analyses take no account of viscous layers at the top and 
bottom of the flow. (The top will have to be capped because 
of large radial pressure gradients.) The main effect of these 
viscous boundary layers is to drive Ekman circulation, which 
flows more rapidly against a weak angular-momentum gradi- 
ent than against uniform rotation. The thickness of the Ek- 



man layer &e 



'v/Q, is small (~ 10 h at point C), 



and the Ekman circulation time C^^^hj ^ftAl ~ 2 s is much 
longer than a typical MRI growth time, so we do not expect 
these layers to be important for stability. To further mini- 
mize their effect, if necessary, one could increase e, use dif- 
ferentially rotating rings at the vertical boundaries, or mod- 
ify the boundary layer by localized Lorentz forces ( [Brahme 
19701). 



A third issue is flnite-amplitude or nonlinear hydrody- 
namical instability in Rayleigh- stable regimes. Few theoreti- 



cal studie s on this subject exist (Serrin 1959; Joseph fc Mun 



Bon 197C). It has been argued from experim ents that a rapid 
Coue tte flow can be nonlinearly unstable (Richard & Zahn 
1999). However, there are indications that such instabilities 



are caused by wall surface defects ( Schultz-Grunow 1959), 
which can be minimized. In fact, it has been shown numer- 
ically that a positive angula r-momentum gradient strongly 
resists nonlinear instability ([Balbus, Hawley, fc Stone 199( ; 



Hawley, Balbus, fc Winters 1999). The outcome depends, 
however, on the amplitude of the initial perturbation and 
the strength of the angular-momentum gradient. These ques- 
tions could be addressed empirically and relatively inexpen- 
sively in prototype experiments using water. 

One would like to predict the nonlinear phases of MRI in 
the laboratory system. Ultimately, nonlinear MRI is closely 
related to other important physics of accretion disks involv- 
ing magnetic fleld, i.e., dynamo processes and jet formation. 



Ongoing three-dimen sional MHD simulations (Kageyama, 



Ji, fc Goodman 2001) will provide useful insights here. In- 



deed, an important purpose of the experiment is to be a 
testbed for MHD codes since, lacking detailed observational 
constraints, theorists depend upon computer simulations to 
understand MRI-driven turbulence. 



fc Miyama 199^ [Fleming, Stone, fc Hawley 2000[ ). On the 



other hand, the innermost (and therefore most energetic) 
parts of accretion disks often encounter a vertical fleld due 
to the central compact object. The works cited above flnd 
that in the presence of an imposed vertical fleld, turbulence 
is driven by axisymmetric modes and persists to higher re- 
sistivity, probably into the experimentally accessible regime. 

In summary, we have used linear stability analyses to 
explore the prospects for magnetorotational instability in a 
magnetized Couette flow. We flnd that MRI can be achieved 
in a moderately rapidly rotating table-top apparatus using 
an easy-to-handle liquid metal such as gallium. Auxiliary 
experiments with an inexpensive nonmagnetic fluid, such as 
water, will be valuable both as prototypes and as controls to 
distinguish MRI from nonlinear hydrodynamic instabilities. 
Onset and dynamics of MRI can be detected by torque mea- 
surements of cylinders and magnetic sensors placed around 
the annulus. Ultrasonic imaging may also be possible. If suc- 
cessful, this will be a rare example of an astrophysical process 
that can be studied in the laboratory. 
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